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ABSTRACT 

The  generation  and  evolution  of  small  amplitude  long  wavelength  travelling  distur¬ 
bances  in  rotating  disk  flow  is  the  subject  of  this  paper.  The  steady  rotational  speed  of  the 
disk  is  perturbed  so  as  to  introduce  high  frequency  oscillations  in  the  flow  field.  Secondly, 
we  introduce  surface  imperfections  on  the  disk  such  as  roughness  elements.  The  interac¬ 
tion  of  these  two  disturbances  will  generate  the  instability  waves  whose  evolution  is 
governed  by  parabolic  partial  differential  equations  which  are  solved  numerically.  It  is 
found  that,  for  the  class  of  disturbances  considered  here  (wavelength  on  the  order  of 
Reynolds  number),  eigensolutions  exist  which  decay  or  grow  algebraically  in  the  radial 
direction.  However,  these  solutions  grow  only  for  frequencies  larger  than  4.58  times  the 
steady  rotational  speed  of  the  disk.  The  computed  receptivity  coefficient  shows  that  there  is 
an  optimum  size  of  roughness  for  which  these  modes  are  excited  the  most.  The  width  of 
these  roughness  elements  in  the  radial  direction  is  about  .lr0*  where  rQ*  is  the  radial 
location  of  the  roughness.  It  is  also  found  that  the  receptivity  coefficient  is  larger  for  a 
negative  spanwise  wavenumber  than  for  a  positive  one.  Typical  wave  angles  found  for 
these  disturbances  are  about  -26®. 
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1.  INTRODUCTION 

In  this  paper,  we  investigate  the  receptivity  and  non-parallel  stability  problem  for 
travelling  disturbances  in  a  rotating  disk  flow.  A  significant  feature  of  the  boundary- 
layer  flow  on  the  rotating  disk  is  that  its  velocity  component  in  the  radial  direction,  the 
crossflow  velocity,  has  an  inflectional  point.  This  causes  an  instability,  often  called 
crossflow  instability,  first  noticed  by  Smith  (1946)  and  investigated  theoretically  and 
experimentally  by  Gregory  et  al.  (1955).  When  the  crossflow  component  is  combined  in  a 
particular  direction  with  the  velocity  component  in  the  azimuthal  direction,  they  form  a 
mean  velocity  profile  which  has  an  inflection  point  at  which  the  velocity  is  zero.  This 
permits  an  inviscid  neutral  disturbance  with  zero  frequency  which  connects  to  unstable 
modes  at  other  wavenumbers  (see  Stuart  in  Gregory,  Stuart  &  Walker  (1955)).  The 
stationary  disturbances  appear  as  vortices,  called  crossflow  vortices,  which  rotate  in  the 
same  direction.  This  phenomena  is  also  observed  in  other  flow  geometries  such  as 
rotating  cones  and  swept  wings,  and  it  is  this  latter  practical  application  which  generates 
the  interest  in  the  problem  since  it  is  found  that  the  crossflow  instability  dominates  the 
boundary-layer  transition  process  near  the  leading  edge  of  a  swept  wing. 

In  the  rotating  disk  flow,  the  crossflow  vortices  spiral  outward  at  an  angle  whose  mag¬ 
nitude  is  approximately  10°-14°,  with  respect  to  the  azimuthal  direction  (Fig.  1),  and  these 
vortices  first  start  to  appear  at  a  Reynolds  number  of  about  286.  Wilkinson  and  Malik 
(1983),  using  hot-wire  techniques,  mapped  out  the  complete  wave  pattern  on  the  disk  and 
found  that  the  stationary  disturbances  originate  from  isolated  roughness  sites  on  the  disk. 
Mack  (1985)  computed  the  wave  pattern  observed  in  the  Wilkinson-Malik  experiment 
using  the  linearized  stability  equations  and  assuming  a  white  spectrum  at  the  source  of  the 
wave  pattern.  Malik  (1986)  calculated  the  neutral  curve  for  stationary  disturbances,  and 
found  numerically  another  stationary  viscous  mode  which  corresponded  to  zero  wall- 
shear  stress  of  the  mean  velocity  profiles.  At  high  Reynolds  numbers  on  the  lower  branch 
of  the  neutral  curve,  the  wavenumber  a  behaves  like  a  -  const/R1/2,  and  at  the  upper 
inviscid  branch  they  behave  like  a  ~  const.  Here  a  is  the  nondimensional  wavenumber  in 
the  radial  direction,  which  is  nondimensionalized  by  the  length  scale  (v/D0)1/2,  and  R  is 
the  Reynolds  number,  defined  by  R  =  r,‘(Do  /  v)1/a,  where  v  is  the  kinematic  viscosity,  Q0\s 
the  angular  velocity  of  the  disk,  and  r(*  is  the  radius  where  the  local  stability  analysis  was 
done. 

Hall  (1986)  developed  an  asymptotic  analysis  to  investigate  the  two  branches  of  the 
neutral  curve.  He  showed  that  the  inviscid  mode  is  characterized  by  a  two-layer  structure. 
The  upper  deck  is  governed  by  the  Rayleigh  equation,  and  a  thin  layer  near  the  wall  is 
needed  to  satisfy  the  no-slip  boundary  condition  at  the  wall.  The  lower  branch  neutral 
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curve  is  characterized  by  a  triple-layer  structure  and  has  a  long  wave  length.  This  mode 
has  its  structure  fixed  by  a  balance  between  viscous  and  Coriolis  forces.  At  high  Reynolds 
numbers,  the  lower  branch  neutral  waves  are  inclined  at  about  39°  to  the  radial  direction. 
It  is  also  observed  that  the  formal  asymptotic  theory  results,  and  the  linear  stability 
calculations  based  on  parallel-flow  assumption,  agree  for  a  large  range  of  finite  Reynolds 
number. 

In  a  recent  paper  by  Balakumar  &  Malik  (1990)  (referred  to  as  BM  below),  linear 
stability  characteristics  of  stationary  and  travelling  disturbances  in  a  rotating  disk  flow 
were  presented.  Figure  (2)  shows  the  neutral  stability  curves  in  the  radial  wavenumber  a 
vs  Reynolds  number  R  plane  for  different  nondimensional  frequencies.  Here  the  nondi- 
mensional  frequency  fit)  is  defined  as  the  ratio  of  the  dimensional  frequency  of  the  distur¬ 
bance  of  to  the  rotational  speed  of  the  disk  £2q.  For  a  detailed  discussion  of  this,  the  reader 
is  referred  to  BM.  In  this  paper,  we  are  concerned  with  the  lower  branch  neutral  curve  at 
high  frequencies.  It  is  seen  from  Fig.  (2)  that  at  high  frequencies,  ct)  =  6.0  and  7.9,  the 
wavenumber  a  varies  like  R1.  The  waves  for  these  modes  are  inclined  at  negative 
angles,  see  Fig.  (1).  For  (0  =  7.9,  this  angle  is  about  -35.32°. 

The  stability  results  in  BM  have  been  obtained  from  linear  stability  theory  with  paral¬ 
lel  flow  assumptions  being  used.  In  the  rotating  disk  case,  this  assumption  amounts  to 
locally  replacing  the  variable  radius  r*  by  the  constant  radius  r*.  This  approach  is  justi¬ 
fied  if  the  wavelength  is  much  smaller  than  the  radius  r*.  This  requirement  is  violated 
when  the  wavenumber  is  of  the  order  of  R'1,  and  the  validity  of  the  Orr-Sommerfeld  type 
approach  becomes  questionable.  This  is  similar  to  the  situation  which  exists  in  the  Gortler 
problem,  Hall  (1982,1983).  It  was  shown  by  Hall  that,  at  0(1)  wavenumbers,  the  Gtirtler 
instability  is  governed  by  a  parabolic  type  partial  differential  system.  In  this  paper  we 
investigate  this  R1  BM  mode  of  instability  at  high  frequencies  and  allow  for  nonparallel 
effects  in  a  self-consistent  manner. 

We  formulate  the  problem  as  a  receptivity  problem  similar  to  that  analyzed  for  the 
Gortler  case  by  Denier,  Hall,  &  Seddougui  (1990).  We  perturb  the  steady  rotational  speed  of 
the  disk  so  as  to  introduce  high-frequency  oscillations  in  the  flow  field.  Secondly,  we 
introduce  imperfections  on  the  surface  as  roughness  elements.  The  interaction  of  these  two 
flow  effects,  the  unsteady  one  due  to  the  oscillations  of  the  disk  and  that  introduced  by  the 
roughness,  will  generate  the  instability  waves  whose  evolution  is  governed  by  parabolic 
partial  differential  equations.  These  partial  differential  equations  are  solved  numeri¬ 
cally.  In  Section  2  we  describe  the  formulation  of  the  receptivity  problem  for  the  rotating 
disk  flow,  the  results  are  discussed  in  Section  3,  and  the  conclusions  are  presented  in 
Section  4. 
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Consider  an  infinite  disk  rotating  about  its  axis  with  angular  velocity  Qq  (Fig.  (1)). 
We  take  cylindrical  coordinates  r*,  6,  z  with  z  =  0  being  the  plane  of  the  disk  and  assume 
that  the  fluid  occupies  the  half  space  z*  >  0.  Let  p,  u,  v,  w  denote  the  steady-state  values  of 
pressure  and  velocity  in  the  r* ,  0,  z  directions,  respectively,  in  the  rotating  coordinate 
frame.  Von  Karman's  exact  solution  of  the  Navier-Stokes  equations  for  steady  laminar 
rotating-disk  flow  takes  the  form 

u  =  r*£20F(z )  ,  v=r'n0G(z)  ,  w  -{vQ^)1'2 H(z)  ,  p  =  pv£20niz )  (2.1) 

where  z  =  z*(fl0/  v)1/2.  The  Navier-Stokes  equations  reduce  to  the  following  equations  for 
F,  G,  H  and  P: 

F2-(G+1)2+F'H-F"  =  0  (2.2) 

2F(G  +  1)  +  G'H -G"  =  0  (2.3) 

n'+HH’-H*  =  0  (2.4) 

2F  +  H'  =  0  (2.5) 

where  the  prime  denotes  differentiation  with  respect  to  z.  The  boundary  conditions  are: 


F=0,  <3  =  0,  *7=0 


(z  =  0) 


F=  0  G=- 1 


(z  — »  «>) 


In  the  following  analysis  the  variables  are  nondimensionalized  by: 
velocity  -  , 


length 


pressure  -  pr? Go  » 


(  V/2  i 

1  —  1  1 

V  A>  J  Gq 


The  Reynolds  number  is  defined  by 
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where  r*  is  a  typical  radial  length  scale. 

Let  us  write  the  instantaneous  nondimensional  velocities  u,  v,  u>  and  pressure  p  as 
u(r,9,z,t)  =  -^F(z)  +  U(r,9,z,t)  ,  ( 

it 


v(r,  d,z,t)  =  -^G(z)  +  V(r,9,z,t)  , 

il 


(2.10) 


w(r,9,z,t)  =  —  lU.z)  +  W(r,f',z,t)  , 
R 


(2.11) 


p(r,9,z,t)  =  —^n(z)  +  P{r,9,z,t)  , 


(2.12) 


where  U,  V,  W  are  nondimensional  perturbations  from  the  steady  mean  flow.  We  substi¬ 
tute  these  expressions  into  the  Navier-Stokes  equations  in  the  rotating  frame  and  obtain  the 
following  nonlinear  perturbation  equations  for  U,  V,  W,  and  P. 


OU  r  9U  GdU  HdU  F rr  r  2 V2  T1dU  V  OU  wdU 

+  -F-— — f/  +—WF-—(G+  1)V - +  U-—  +  — —  +  W-— 

dtR  drR  dd  RdzRR  R  r  dr  r  dd  dz 


dP  1  \d2U  1  d2U  1 
"  dr  +  R\dr2+r2  d02+r 


du  dHj_  u_  2  jv 

dr  +  dz2  r 2  r2  do  ' 


(2.13) 


dV  r  dV  GdV  H  dV  F..rwn,  2  UV  T1dV  V  dV  mdV 

— -  +—F-—  +—-— — V  +  —WG+ —{G  +  DU  + +  17— -  + —  +  W— 

dt  R  Or  R  09  RdzRR  R  r  dr  r  09  dz 


1  dP_  1  d2V  1  d2V  IdV  £V__V_  2  OU 

r  09  +  R  dr2  +  r2  092  +  r  dr  +  dz2  r2  +  r2  99 


(2.14) 


dW  r  „dW  GdW  H  OW  1  rrdW  V  OW  ,,.dW 
dt  R  dr  R  09  R  dz  R  dr  r  09  dz 

dP  1  f(?2W  1  d2W  1  OW  <92W] 

dz  R  dr2  r 2  dO2  r  dr  dz2  J 


(2.15) 


u  [  ou  |  i  ov  t  OW 

r  dr  r  09  dz 


(2.16) 


Here  r  is  the  nondimensional  radius  and  t  is  the  nondimensional  time.  We  must  solve 
(2.13)  -  (2.16)  subject  to  the  no-slip  condition  at  the  wall,  whereas  sufficiently  far  away 
from  the  wall  we  insist  that  the  disturbances  vanish.  We  are  going  to  investigate  the  high 
frequency  lower-branch  instability  mode  for  which  the  wavenumber  varies  like  R  l  at  high 
Reynolds  numbers  as  discussed  in  BM. 
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We  introduce  a  disturbance  into  the  steady  motion  first  by  perturbing  the  speed  of 


rotation  of  the  disk  to 
fl  =  flb(l+4e-ia,v)  , 


(2.17) 


and  second  by  introducing  a  hump  on  the  surface  of  the  disk  (Fig.  1).  Here  o>*  is  the  dimen¬ 
sional  frequency  of  the  oscillations  which  may,  for  example,  simulate  an  acoustic  distur¬ 
bance.  The  shape  of  the  hump  is  defined  by 


(  \1/2  (  *  \ 
'■■(i)  "(H  ■ 


(2.18) 


where  /  is  a  function  of  (r*  /  r* )  and  0;  8  and  A  are  small  scalar  parameters  with  8,  A  «  1. 
In  nondimensional  form  (2.17)  and  (2.18)  become 


a=n0(i+Ae~i<o,Rt)  , 


(2.19) 


(2.20) 


where  <o  is  the  nondimensional  frequency  defined  by  CD  =  (o*  /  Q0.  Since  we  are  seeking 
solutions  which  are  linear  in  8,  we  can  look  for  Fourier  series  solutions  in  the  6  direction. 
We  expand  fir!  R,8)  in  Fourier  series  as 


(2.22) 


This  will  permit  us  to  perform  the  analysis  in  the  Fourier  space  fl.  We  will  assume  that 
Reynolds  number  R  is  large  and  seek  solutions  which  are  linear  in  8  and  A.  The  perturba¬ 
tions  U,  V,  W  and  P  are  written  as 

U  =  A  ^-F(z)e~Uo,R‘  ^,z\em 

R  R  \R  ) 

+A5u^,x 

+ higher  order  terms  (2.23) 

V  =  AjjG(z)e-ie,/Rt 

+  higher  order  terms  >  (2.24) 
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W  =  A^H(2)e-ia,/R‘ 

R  R  \R  J 

+A.5±u>[j,zy-ia,'R,eif,B 

+  higher  order  terms 

We  identify  ( F,G,H,P )  as  the  flowfield  induced  by  the  oscillations  of  the  disk,  ( F,G,H,P ) 
as  the  flowfield  induced  by  the  hump,  and  ( u,v,w,p )  as  the  flowfield  produced  by  the 
interaction  between  the  first  two  flow  fields.  Let  us  assume  that  the  factor  r  -r  I R  is  0(1) . 
Substituting  the  above  expressions  (2.23)  -  (2.26)  into  the  equations  (2.13)  -  (2.16)  and 
collecting  terms  of  the  same  order  we  obtain  a  sequence  of  equations  for  the  induced  fields. 


At  order  A,  we  obtain 

(-»©  +  2  F)F  +  HF’ -2(G  +  1)G  +  F3f -F*  =  0  (2.27) 

{-ia  +  2F)G  +  HG'  +  2  (G  +  1)F  +  G'H  -  G"  =  0  (2.28) 

{-ia>+H')H+HH'+P-H~  =  0  (2.29) 

2F  +  H’  =  0.  (2.30) 

The  boundary  conditions  at  the  solid  boundary  are 

F(0)  =  hT(0)  =  0  (2.31) 

G(0)  =  1  .  (2.32) 


Far  away  from  the  surface  of  the  disk 
gives 

F'-fjF  =  0 

&-pG  =  0 

where 

H(oo)-Vh2(oq)-4zq? 

M  2 


the  disturbance  should  decay  exponentially.  This 

(2.33) 

(2.34) 

(2.35) 


(2.25) 


(2.26) 
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and  the  real  part  of  fi  is  less  than  zero.  For  order  5  we  obtain  the  following  equations  for 
F,G,H  and  P: 


(2F  +  iPG)F  +  rF^~  +  H~--2(G+l)G  +F'H  =  , 

dr  dz  dz * 

(2.36) 

(2F  +  ipG)G  +  rF^  +  H^  +  2(G  +  l)F  +  G'H  =  ^  , 

dr  dz  dz 1 

(2.37) 

rF~-+ipGH  +  H^-  +  H'H  =  -^-+£K 
dr  dz  dz  dz2  * 

(2.38) 

2F  +  r~  +  ipG+~  =  0  . 
dr  dz 

(2.39) 

The  appropriate  boundary  conditions  are 

F  (0)  =  -F’(0)f  (r,p) 

(2.40) 

G(0)  =  -G'(0)f(r,P) 

(2.41) 

H(  0)  =  0 

(2.42) 

The  solution  of  Eq.  (2.36-2.42)  is: 

F(z)  =  -F'(z)f(r,P)  , 

(2.43) 

G(z)  =  -G'(z)f(,r,P)  , 

(2.44) 

H(z)  =  ■  2f(r,p)  +  r^F(z)  +  ipG(z)hr,p)  , 

(2.45) 

and  P(z)  can  be  obtained  from  Eq.  (2.38).  At  order  AS,  we  obtain  the  equations  for  the 
induced  disturbances  u,v,w and  pi 

~  -  (-i<a  +  ipG  +  F)u-rF^-H~- FwF'  +  2  (G  +  1)0  =  |-2rGG  +  FF  j  F  +  r  ^  j 

+FFF  +  ipGFr  +  rH~  +  rH^~\  ,  (2  46) 

dz  dz  j 

^-(-ia>  +  ipG  +  F)v-rF^-H~-rG'w-2(G+V>a  =  rl*'G+FG)  +  rF(G+r  — 

dz  dr  dz  '  '  \  or 
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(2.47) 


+rFG  +  ipGGr  +  rH^  +  rH > 

OZ  oZ 

dp  d2w  ,  .  .or%  .  __  dw  TT  dw  dH  -oATT  rrdH  77  dH  ao. 

~= — -T  +  (-m  +  i{iG  +  H')w  +  rF-—  +  H—  +  rF—  +  if)GG  +  H—  +  H—  >  (2.48) 
oz  oz  of  oz  of  oz  oz 


u  +  r  —  +  t/fo  +  r  —  =  0  . 

of  oz 


(2.49) 


The  appropriate  boundary  conditions  are 


u(0)  =  ~F'(0)rf(r,p )  , 


(2.50) 


m=-G’mhr,p)  , 


(2.51) 


iKO)  =  0 


(2.52) 


It  is  easy  to  verify  that  if  we  write 


u  =  -F'(z)rf(?,P)  +  ul  , 


(2.53) 


v  =  -G\z)rf(r,p)  +  vl  , 


(2.54) 


u)=  2f  +  r^^F(z)  +  ipGf  +  wx  , 


(2.55) 


then  Eq.  (2.46),  (2.47)  and  (2.49)  take  the  form 


■^-^1-  =  (-i<a  +  i)3G  +  F)u1  +  rF-^-  +  fl'-^-  +  rF,u;i -2(G  +  Dt^  , 
oz*  or  oz 


(2.56) 


^  =  ( -i(o  +  ipG  +  F)u,  +  rF  %  +  H  % + fG’wx  +  2(G  + D«i  , 

OZ*  of  oz 


(2.57) 


_  dw,  _<?w,  ,a  n 

r Ik  +  ui  +  r#+lPvi  =  °  > 


(2.58) 


with  the  boundary  conditions 
u1(r,0)  =  0  , 


(2.59) 


i>i(r,0)  =  0  , 


(2.60) 


wx(r,Q)  =  -ipf{r,P)  . 


(2.61) 


The  boundary  conditions  at  the  far  field  are 


v[  -jhjl=Q 


(2.63) 


where 


.  H(o°)-Jmoo)2-i((0+p) 

2 


(2.64) 


and  the  real  part  of  fi  is  less  than  zero.  Equation  (2.48)  can  be  integrated  to  solve  for 
pressure  p.  We  observe  that  the  above  system  of  equations  is  of  the  parabolic  type.  The 
inhomogenity  appears  as  the  normal  boundary  condition  u;,(r,0)  =  Therefore,  if 

we  specify  <o,  p  and  f(r,P)  we  can  integrate  these  equations  marching  in  r .  The  equations 
were  solved  using  the  two-point  fourth-order  compact  scheme  in  the  z  direction  and  second- 
order  upwind  scheme  in  r  direction  (Spall  &  Malik  (1989)). 

We  note  that  the  system  (2.56-2.58)  permits  eigensolutions  of  the  form 


(2.65) 


where  yis  a  complex  number  and  (uo,v0,u>0)  satisfy, 

^  =  (-i(o  +  ipG  +  F)u0  +  iFu0+H^-  +  w0F/-2(G+l)v0  ,  (2.66) 

dz  dz 

»2  i 

— V-§-  =  (-ioi  +  ipG  +  F)v0  +  yFvQ  +  H— —+w0G'+2(G+l)ui}  ,  (2.67) 

dz  dz 

^7^  +  «o  +  yu o  +  iPv o  =  0  ,  (2.68) 


with  the  boundary  conditions  Uo(0)-  t>o(0)  =  u;o(0)  =  0  and  u,  v  and  w  decay  at  infinity. 
Therefore,  yis  the  eigenvalue  of  the  system  and  determines  whether  the  disturbance  grows 
or  decays  algebraically  in  the  radial  direction.  If  in  the  parallel  flow  linearized  stability 
equations  (2. 16-2. 19)  of  BM,  we  replace  a  =  a0  /  R,  P  =  P0/  R,  at  =  a>  /  R  and  neglect  all  the 
l/f22  terms,  we  obtain  the  same  equations  as  above  (2.66-2.68)  with  ia  replaced  by  y. 
Therefore,  the  solutions  obtained  for  the  high  frequency  (UR)  mode  in  the  earlier  paper  are 
to  be  interpreted  as  algebraic  and  not  as  exponentially  growing/decaying  solutions.  In 
fact,  the  linearized  Navier-Stokes  equations  in  a  rotating  coordinate  system  (equations 
(2.11-2.14)  in  BM)  have  solutions  of  the  form 
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g.rf{*w+J^+-.+4£>+..|  , 


r  s  ,  s  V,(z)  .  Vn(z ) 

i>  =  ryji>o(z)  +  -L-+»-+-J^r  +  -  -  , 


-y[iZ'0(z)  ,  iii(z)  ,  ^»(z)  ,  1 
~  I  r  r2  r2n+1  |  ’ 


.  rr  \f>0  (z)  .  ji,(z)  ,  .  £„(z)  .  1 

p—r— +_ !^-+-+— ^ *"•  ' 


This  solution  exists  for  all  the  frequencies  co  and  all  azimuthal  wavenumbers  /?.  However, 
we  will  see  in  the  next  section  that  the  real  part  of  yis  positive  only  for  frequencies  o>  >  4.58 
and  for  a  certain  range  of  /}  values,  and  thus  the  algebraic  instability  is  not  present  for 
lower  frequencies. 


3.  RESULTS 

First  we  will  present  the  values  of  y obtained  from  the  eigenvalue  system  (2.66-2.68). 
The  results  are  depicted  in  Table  1. 

Table  1 


+4.0 

+0.0 

(-.50388,15.318) 

-1.0 

(-.48892,13.176) 

-2.0 

(-.49032,10.829) 

-3.0 

(-.54533,8.089) 

+4.6 

-2.0 

(-.03848,13.484) 

-1.0 

(-.00918,15.783) 

0.0 

(.00422,17.916) 

1.0 

(.00921,19.936) 

2.0 

(.00753,21.874) 

3.0 

(.00186,23.748) 

4.0 

(-.00675,25.570) 

+7.9 

-8.0 

(-1.3621,9.593) 

-7.0 

(1.1622,14.652) 

-6.0 

(1.6793,18.261) 

-5.0 

(1.9842,21.249) 

-4.0 

(2.2082,23.919) 

-2.0 

(2.5397,28.716) 

0.0 

(2.7884,33.074) 

2.0 

(2.9893,37.156) 

4.0 

(3.1582,41.045) 
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We  observe  that  the  real  part  of  yis  positive  for  frequencies  larger  than  4.58,  however,  this 
occurs  only  for  a  certain  range  of  P  values.  For  co  =  4.6,  real  (y)  is  positive  in  the  range  of  /? 
=  0.0  to  3.0  and  for  co  =  7.9,  this  range  is  from  -7.0  to  40.  These  results  are  in  agreement  with 
the  BM  results  reproduced  in  Figure  2  which  shows  that  UR  mode  instability  is  present  only 
for  frequencies  higher  than  4.58.  The  unstable  modes  computed  by  BM  for  frequencies 
lower  than  4.58  do  not  belong  to  this  class  of  instability. 

Now  we  present  the  results  for  the  receptivity  coefficient  associated  with  the  UR  mode 
discussed  above.  First  we  consider  the  case  where  f(.r,fi)  takes  the  form 

f(r,p)  =  e~ *  , 

with 

x  =  Vcr(r-r0)  . 

Here  r0  is  the  location  of  the  hump  and  will  be  taken  as  r0  =1,  which  is  equivalent  to  defin¬ 
ing  the  r*  based  on  the  radial  location  of  the  hump.  The  parameter  a  determines  the  spread 
of  the  hump  in  the  radial  direction. 

Figures  3,  4  and  5  show  the  downstream  development  of  the  amplitude  of  ux,  vl  and  wx 
(Eqs.  (2.56-2.58))  for  the  case  (o  =  7.9,  P  =  -7  and  a  =  40.0.  In  these  figures,  the  vertical  axis 
shows  the  normal  coordinate  z  and  the  different  curves  are  plotted  at  different  r  =  .6,  .7,  .8, 
.9,  1.0,  1.1,  1.25,  1.3,  1.4,  1.6  and  1.9.  We  notice  that  the  initial  form  of  the  radial  distur¬ 
bance  velocity  component  ux  shows  two  maxima,  and  the  lower  maximum  disappears 
while  the  amplitude  of  the  upper  maximum  increases  continuously  with  r.  The 
distribution  of  the  amplitude  of  the  normal  velocity  component  wx  shows  that  at  the  edge  of 
the  boundary  layer  it  approaches  a  constant  value  which  increases  with  increasing  r . 

In  Figs.  6,  7,  8,  and  9  we  plot  the  maximum  amplitude  of  the  velocity  components  ux,  vx, 
and  wx  against  the  coordinate  x  for  various  values  of  a  with  p  =  -7  and  (o  -  7.9.  For  ux  we 
take  the  maximum  amplitude  at  the  outer  maximum  point  and  for  wx  we  take  the  maxi¬ 
mum  amplitude  as  the  amplitude  at  the  edge  of  the  boundary  layer.  We  observe  two  differ¬ 
ent  types  of  curves  for  a  fixed  a.  For  an  x  less  than  about  .5,  the  amplitude  pattern  appears  to 
be  Gaussian  shaped,  and  beyond  that,  the  amplitude  increases  with  increasing  x.  We 
identify  that  the  first  part  is  due  to  the  direct  influence  of  the  wall  and  the  second  growing 
part  is  due  to  the  instability.  Figure  9  is  the  same  as  Fig.  7  and  the  curve  A  =  Cx(p,o,(o)rr  is 
plotted  for  a=  10,  20,  100,  and  200.  Here  y  is  the  eigenvalue  obtained  from  Eq.  (2.66-2.68)  for 
0)  =  7.9  and  P  =  -7  and  is  equal  to  (1.1622,  14.652).  Cx{p,o,(o)  is  a  measure  of  the  receptivity 
coefficient  of  these  modes  to  the  excitation  we  considered  in  this  paper.  From  Fig.  9  it  is 
clear  that  there  exists  downstream  an  algebraic  instability  of  the  form  rf  and  there  are  no 
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other  eigenfunctions  which  exist  for  the  system  (2.56-2.58).  The  receptivity  coefficient 
C,(/J,<7,fl>)  is  given  in  Table  2  for  (o  =  7.9  and  P  =  -7  for  different  a  values. 


Table  2 


<T 

Ci(0=-7,<T,fi)=7.9) 

1C,  1 

5 

(  .00135, .00324) 

.00352 

10 

(-.00391, .00983) 

.01058 

20 

(-.04381..01251) 

.04556 

40 

(-.11258, .01124) 

.11314 

100 

(-.16154, .02733) 

.16384 

200 

(-.14827, .03586) 

.15254 

1000 

(-.08 145, .02089) 

.08409 

2000 

(-.05812,-01693) 

.06054 

The  receptivity  coefficient  increases  with  <x  up  to  a  -  100  and  then  decreases.  Therefore,  we 
conclude  that  there  exists  an  optimum  size  for  which  instabilities  are  excited  the  most.  In 
Fig.  10,  we  plot  the  maximum  amplitude  of  vx  for  (0  =  7.9,  a  =  100  for  different  P  values.  The 
figure  shows  that  the  waves  are  excited  for  the  range  of  P  values  P  =  -7  to  p  =  20.  The 
receptivity  coefficients  are  larger  for  negative  P  values  than  for  positive  P  values. 

Next,  we  will  present  the  cumulative  wave  pattern  produced  downstream  of  the  rough¬ 
ness  site.  We  will  consider  the  case  0)=  7.9  and  o-  100  and  the  roughness  element  defined 
by 

f(r,6)  =  e-x2S(0 )  , 

where 

x  =  4o{r  -  L0)  , 

and  S(0)  defines  the  shape  in  the  azimuthal  direction.  We  consider  the  triangular  shape 
defined  by 

S(6 )  =  — (0  +  a)  ;  -oie<,  0 
a 

=  -£■(6-0)  ;  OS0S6 

0 

where  a,  b  are  positive  numbers  which  are  in  the  range  0  <  a,  b  <  n.  If  a  equals  b  we  get  the 
symmetric  shape;  otherwise,  we  get  the  asymmetric  case.  We  will  present  the  results  for 
the  cases  a  =  b  =  .05.  The  Fourier  transform  of  S(0)  is 
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S(0)  =  £S(/ 3)eip9 


where 


S(/3)  =  f-(a+6)  ;  0  =  0 

Zit 


1  1  fl-e,/fa  1-e"1^ 


0*0  . 


Figure  1 1  shows  the  contour  plot  in  (r ,  0)  plane  of  the  real  part  of  t>4 ( r,9,z,t )  at  the  height  z  = 
1.5  for  the  case  to  =  7.9  and  a  -  100,  a  -  .05  and  b  =  .05.  The  roughness  element  is  located  at 
ra  =  1  and  is  denoted  by  the  point*.  Near  the  roughness  we  see  the  effect  of  the  roughness, 
and  downstream  of  the  roughness  we  observe  the  instability  waves  emerging  from  the 
wake.  Contours  of  zero  amplitude  are  highlighted  by  thick  lines.  These  waves  are 
inclined  at  angles  of  about  -26°.  The  wavelength  in  the  radial  direction  is  about  1.0  and  the 
number  of  waves  in  the  azimuthal  direction  0  is  about  6. 

In  Fig.  12  we  plot  the  maximum  amplitude  of  IV0  for  to  =  4.6  for  different  0  and  cr 
values.  We  see  that  there  is  only  weak  instability  for  this  frequency  for  certain  values  of  0. 
Figures  13  and  14  show  the  results  for  frequencies  (O  -  4.0  and  0.0.  For  these  frequencies, 
all  the  disturbances  decay  downstream.  All  these  conclusions  agree  with  our  previous 
calculations  for  yand  the  results  of  BM  that  the  1 JR  instability  is  present  only  for  a>  >  4.58. 


4.  CONCLUDING  REMARKS 

The  receptivity  and  non-parallel  stability  of  travelling  disturbances  of  long 
wavelength  (on  the  order  of  Reynolds  number  R)  are  investigated.  The  non-parallel 
effects  are  included  in  a  self-consistent  manner.  It  is  shown  that  eigensolutions  exist 
which  grow  or  decay  algebraically  in  the  radial  direction.  It  is  established  that  the 
solutions  obtained  for  the  high  frequency  (1/R)  mode  in  the  earlier  paper  of  BM  are  to  be 
interpreted  as  algebraic  and  not  as  exponentially  growing/decaying  solutions.  These 
algebraic-type  solutions  exist  for  all  frequencies  (0  and  all  spanwise  wavenumbers  0. 
However,  growing  disturbances  exist  only  for  frequencies  larger  than  4.58  times  the  disk 
rotational  frequency  and  for  a  certain  range  of  0  values  which  is  in  agreement  with  the 
results  of  BM. 

The  receptivity  calculations  show  that  there  exists  an  optimum  size  of  roughness  for 
which  instabilities  are  excited  the  most.  For  example,  the  receptivity  coefficient  C\  for 
frequency  (o  =  7.9  with  spanwise  wavenumber  fi  =  -7  is  .01,  .16  and  .084  when  a-  10, 100  and 
1000,  respectively.  Gaussian  shaped  roughness  elements  with  a  =  100  correspond  to  a 
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roughness  element  with  a  spread  of  about  .lrg  in  the  radial  direction,  where  r0*  is  the 
location  of  the  roughness.  This  shows  that  very  narrow  or  broad  roughness  elements  will 
not  excite  these  modes  efficiently.  It  is  also  observed  that  the  receptivity  coefficients  are 
larger  for  negative  values  than  for  positive  fi  values.  Hence,  if  these  modes  are  excited 
they  will  align  along  negative  angles. 

The  cumulative  wave  pattern  produced  from  an  isolated  roughness  element  shows  that 
the  waves  are  inclined  at  about  -26°  and  the  wavelength  in  the  radial  direction  is  about  1.0. 
In  dimensional  variables,  the  wavelength  1.0  equals  to  r0*  where  the  r0*  is  the  location  of  the 
element.  In  some  experiments  performed  by  Faller  and  Kaylor  (1966),  in  Ekman  layer 
and  rotating  disk  boundary  layer,  the  waves  oriented  at  negative  angles  have  been 
observed  at  finite  Reynolds  numbers.  These  waves  were  rapidly  moving  and  had  long 
wavelengths  compared  to  cross-flow  vortices.  Even  though  there  seems  to  be  some 
similarities  between  these  computed  algebraic  modes  and  the  experimental  observations, 
more  measurements  have  to  be  done  to  draw  any  definite  conclusions. 

Another  remaining  question  is  how  the  lower  branch  exponential-type  solutions  join 
with  this  algebraic  mode.  For  example,  consider  the  neutral  stability  curve  for  frequency 
i o  =  6.0  shown  in  Fig.  2.  The  lower  branch  neutral  curve,  which  is  asymptote  to  the  neutral 
curve  for  stationaiy  disturbances,  is  an  exponential-type  solution.  We  see  that  this  curve 
turns  around  at  a  Reynolds  number  of  3xl06  and  becomes  a  1 /R  mode.  However,  we 
discussed  that  this  1/R  mode  is  actually  an  algebraic-type  mode.  Therefore,  why  does  the 
lower  branch  neutral  curve,  which  is  asymptote  to  the  neutral  curve  for  stationary 
disturbances,  turn  around,  and  what  is  its  subsequent  behavior,  are  the  remaining 
questions. 
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